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We calculate the coupling constants of the decay modes 1 h — > piv, fiw, biw, nir, rj'ir, aiir, f\r\ 
within the framework of the light-cone QCD sum rule. Then we calculate the partial width of these 
decay channels, which differ greatly from the existing calculations using phenomenological models. 
For the isovector 1 h state, the dominant decay modes are p7r,/i7r. For its isoscalar partner, its 
dominant decay mode is a\-K. We also discuss the possible search of the 1 h state at BESIII, 
for example through the decay chains J/^p(tp') — > n 1 + 7 or J/1/1 (1//) — > n 1 + p where -K\ can be 
reconstructed through the decay modes 7ri — s> pn — > tt + tt~it or 7Ti — > /i(1285)7r°. Hopefully the 
present work will be helpful to the experimental establishment of the 1 _+ hybrid meson. 

PACS numbers: 12.39.Mk, 12.38.Lg 

Keywords: Hybrid meson, Light-cone QCD sum rule 



I. INTRODUCTION 

Quark model has been proved to be very successful in the classification of hadrons and calculation of hadron spec- 
trum and their other properties. Yet Quantum Chromodynamics (QCD), which is widely accepted as the fundamental 
theory of the strong interaction, does not prohibit the existence of those hadron states which can not be accommo- 
dated in the conventional quark model. These non-conventional hadrons include multi-quark states (qqqq, qqqqq, • ■ ■ ), 
glueballs (gg, ggg, •■■), and hybrids (qqg). Searching for non-conventional hadrons experimentally and studies on 
their properties have attracted much interests over the past few decades. With these efforts one aims to explore the 
nonperturbative aspects of QCD in the low energy sector. 

Some non-conventional hadrons are totally "exotic" , namely their J PC quantum numbers are excluded by the 
conventional quark model. For example, a conventional qq meson possesses parity P = (— 1) L+1 and C-parity 
C = ( — 1) L+S , where L and S are the orbital angular momentum and spin of its componential quark and antiquark, 
respectively. Apparently, J PC = 0~~, H , 1 h , 2 H etc. are impossible for conventional mesons. Hadrons with these 
J PC quantum numbers are exotic states, which are widely studied since they do not mix with conventional hadrons. 

Recently, COMPASS collaboration observed a resonance with exotic quantum numbers J PC = 1 h at (1660 ± 
lOlgJ MeV/c 2 with a width of (269 ± MeV/c 2 Q. In literature, three isovector J PC = 1~+ exotic mesons, 

namely 7Ti(1400), 7Ti(1600), and 7Ti(2015), have been reported. Several groups observed 7Ti(1400) in the r/TT~ system 
was first reported through a combined study of the r/'ir , fiir , and p°ir channels by VES in Ref. [|| 
and later in the biir final state in Ref. E852 collaboration observed 7Ti(1600) by carrying out the partial wave 
analysis in the 7r + 7r~7r~ final state Q and examining the ri'ir~ final state in the reaction Tr~p — > pri'iT~ Q. They also 
found it in the final states /i7r Q and &i7r Q. The Crystal Barrel collaboration analyzed the reaction pp — > W7r + 7r _ 7r° 
and observed the 7Ti(1600) in the bin channel in Ref. @. E852 also reported another 1 ^ meson 7Ti (2015) 0,1- 

The 1 ^ hybrids have been studied in a few different theoretical schemes. The mass of the lowest-lying 1 ^ hybrid 
meson was predicted to be around 1.9 GeV in the flux tube model [l7| . Llanes-Estrada and Cotanch predicted the 
hybrid mass to be above 2.0 GeV utilizing a QCD inspired Coulomb gauge Hamiltonian [l(J. H.-C. Kim and Y. Kim 
investigated the 1 ^ hybrid within the framework of an AdS/QCD model [ll|. Early studies with the leading order 
QCD sum rule estimated the mass of the 1 ^ hybrid to be (1.6 ~ 2.1) GeV [2l|. The inclusion of the radiative 
corrections changed this estimate to 1.26 ±0.15 GeV [2.2|, (1.6 - 1.7) GeV dlLand 1.81(6) GeV [H. The Lattice 
QCD prediction for the mass of 1 ^ falls into a wide range of (1.5 ~ 2.2) GeV [14l ]. 
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The 1 h hybrids were shown to exist as narrow resonant states in the large iV c limit of QCD 12]. The coupling 
of a neutral hybrid {1,3,5...} h to two neutral (hybrid) mesons with the same J PC and J = were studied in 
the large N c limit of QCD in Ref. 13]. Cook and Fiebig presented a decay width calculation on lattice for the 
channel 1 ^ — > a\n in Ref. ,15]. The partial widths of the ground 1 h hybrid to b\n and fin were predicted to be 
400(120) MeV and 90(60) MeV, respectively in Ref. [l6|. The strong decay properties were explored in Ref. [l8[ 
(IKP) within the same framework. Page, Swanson, and Szczepaniak (PSS) extended the original IKP flux tube 
model and studied the strong decays of hybrid mesons with different J PC quantum numbers, including those of the 
1 h hybrids. Burns and Close [20( compared the flux tube model and Lattice QCD for the 5-wave decay of the 1 h 
hybrid and found excellent agreement. 

Some decay modes of the isoscalar and isovector 1 *~ hybrids were studied using the three-point function sum rule 
[25l ] . In Ref. [26| , the decay widths of the 1 v hybrid were calculated using a three-point function sum rule evaluated 
at the symmetric Euclidean point. The mass of the strange hybrid was also studied in this paper within a light quark 
expansion formalism. The partial width of the channel 1 ^ — > pn was predicted to be rather broad by using the 
three-point function at the symmetric point [2H I27I] . Zhu reexamined the decay channels 1 h — > pn, fin by using the 
light-cone QCD sum rules and reduced the partial width of 1 ^ — > pn significantly [Hj]. 

In this work, we employ Light-cone QCD sum rule (LCQSR) 29] to calculate the various coupling constants of the 
decay modes 1 ^ — > pn, f -\n, b-\ n, nn, ri'n, ai n, /] r\. The LCQSR is different from the traditional Shifman-Vainshtein- 
Zakharov (SVZ) sum rule [301 ] . where the short-distance OPE expansion is made. The OPE expansion of the LCQSR 
is performed near the light cone. With the extracted coupling constants, we calculate the partial widths of these 
modes and compare them with the existing results obtained using other approaches. We also suggest the possible 
search for the 1 h hybrid mesons at BESIII. 

The paper is organized as follows. We illustrate the formalism and derive the LCQSR for the tt\ — > pit coupling 
constant in Sec. UU We present the sum rules for the ni — > fit:, bit: and tt\ — ¥ r\n,v[n coupling constants in Sec. 
IIII and ITVl respectively. The decay widths of the isovector and isoscalar 1 h states are presented in Sec. IVl and 
VT1 respectively. We discuss the possible search of these exotic states at BESIII in Sec. IVIII The last section is a 
short summary. The light cone distribution amplitudes of the pion which are employed in the present calculation are 
collected in the appendix. Readers who are not interested in the technical details may skip Sees. ILTllIVI and go to last 
three sections directly. 



II. SUM RULES FOR THE m ->■ pn COUPLING CONSTANT 

We denote the isovector and the isoscalar J PC = 1 h hybrid meson by ni and fri, respectively. The adopted 
interpolating current for ni reads 



1 
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u(x)^-g s G^ u (x)Yu(x) - d(x)^-g s G; v (x)Yd(x) 



(1) 



The overlapping amplitude f 7Tl between the above interpolating current and 7Ti is defined as 

(01^(0)1^^^))=/^^, (2) 

where rj^ is the polarization vector of 7Ti . 

We consider the following correlation functions in our calculation: 

n p (fc 2 ,p 2 ) = % J d*xe ik - x (*(q)\T{JZ(x) j; if (0)}|0) = l e a ^ s k''q 5 G p (k 2 ,p 2 ) + ■■■ , (3) 

where p — k + q. The interpolating current for the p meson is J^{x) = d(x)j p u(x) and we have (0| J^(0)\p(k, A)) = 
fpm,p£* where is the polarization vector of the p meson. The decay amplitude of the channel tt\ — > pn can be 
written as 

Minx -> pn) = e a ^ s e* a V P ^q s g p , (4) 

where k and q are the momentum of the final p and n, respectively. When k 2 ,p 2 <§; 0, G p (k 2 ,p 2 ) can be calculated 
by operator product expansion(OPE) near the light-cone x 2 = 0, with the n light-cone wave functions as its input. 
Furthermore, G p (k 2 ,p 2 ) can be related to g p by the dispersion relation 



G p (k, P )=l d Sl j dS2J -— w -^-— w -— )+ ^ dsi ___ + y o dS2 ____ + ... i (5) 
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where 



p(si,s 2 ) = f P U 1 m p gpS(si - m 2 p )S(s 2 - m\ x ) + 



(6) 








FIG. 1: The Feynman diagrams for H p (k 2 ,p 2 ). 



ISA A'l 

After invoking the double Borel transformation B^ 1 B p2 2 , we extract the double dispersion relation part of Eq. ([S]): 

f P U imp g p e^/ M2+ ^/ M2 + ■■■ 

- iV (A [ r ] A l ? ] vj Ql] - vj Qa] + vi Ql] + V [ ?* ] )mlM 2 - t£t=[M«o) + <M^o)]<S s 2 G 2 ) , 



V2 



where 



it 



M 



36^ 



2 _ AffAfl 



(7) 



(8) 



and x = 1 — x. Hereafter we ignore the factor e ^>uomjM b ecause m l/M 2 < 0.01 in our calculations. The definitions 
of J r [ a *ls are 



/•MO 

jr[»i] = / ,F(ai, ito, uo — ot\) da\ , 
Jo 

/■So 

J 7 !" 2 ! = / J 7 (u ,a 2 ,u - a 2 )da 2 , 
Jo 
/■So /•« 

/ / ^(ai, a 2 , Q-i — a 2 ) da' 2 dai . 
Jo Jo 

/•So /■« 

/ / ^"(ai, a 2 , a 2 — a[) da[da 2 . 
Jo Jo 



Here the Borel transformation is defined as 



(-k 2 ) n+1 ( _d 
dk 



n— >oo 77,! 
2 ^2 



ir 2 ) 



fe 2 =-nM? 



(9) 



(10) 



The quark propagator used in the OPE of G p (k ,p ) is 



^ W ~ 2ttV + 32^ 2 5sG ^x2 (fT ^ + ^ } 12 192 



We present the Feynman diagrams corresponding to the quark-level calculation of TL p (k 2 , p 2 ) in Fig. [TJ 
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The spectral density p(s\, s 2 ) can be derived from the dispersion realtion Eq. ([5]) after two continuous double Borel 
transformation of G p (k 2 ,p 2 ): 



P { Sl , s 2 ) = B- a \ bp 2 b;i s;i G p {k*y) 



(ii) 



According to quark-hadron duality, we can subtract the contribution of the excited states and the continuum from 
Eq. ([7]) and arrive at 



u m 2 p /M 2 +u ml i /M 2 _ 



ds 



2 e 



- Sl ^ e s 2 a 2 g.i B s 2 B g Gp(fc 2 )p2) 



(12) 



where S01 an d S02 are the continuum thresholds of the mass rules of the p meson and the hybrid state 7r 1; respectively. 
The large mass difference between the ix\ hybrid and the p meson inclines us to work at an asymmetric point of the 

Borel parameter M 2 and M|, leading to a sophisticated subtraction of the continuum contribution (34|. The terms 

11 

of B£i B°i G p (k 2 ,p 2 ) have general form cu n (M 2 ) n = co™ j{p\ + er 2 )" l+n . Here we assume m, n > to illustrate the 
procedure of the continuum subtraction. 



dsi I ds 2 e -^i e - S2<J2 B s a \Bo 



J 

Sol /*S02 

dsi I 
Jo 

2 / ds 
Jo 



(<7i + cr 2 ) m +' i 



ds 2 e - Siai e- S2rT2 



r(m + n) 
{*+- 



dS(si - s 2 ) 



9s 1 



y rn\ 8 01 
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2 m r(r7i + n) V9s_ 



m+n— 1 



5(2s_) 



(13) 



where s + = (si + S2)/2, s_ = (s 2 — s i)/2, \jM 2 _ — 1/M 2 — 1/Mf and we assume S01 < S02 ■ /n(^) is the subtraction 
function defined as 



2 

/„(x) = l-e-^^. 



(14) 



i=0 



The 1 + hybrid has not been firmly established experimentally. Both theoretical predictions and experimental 
measurements suggest that the mass of n\ falls within the range 1.6 ~ 2.0 GeV. In this work the mass of m is 



taken to be 



1.6, 1.8 and 2.0 GeV. We adopt f 7Tl = 0.15 GeV in our numerical analysis [231 ] . The n decay 



constant = 131 MeV. The mass and the decay constant of the p meson are m p = 0.77 GeV and f p = 0.216 GeV. 
p„ = m 2 J(m u + m d ) = (1.573 ± 0.174) GeV is given in Ref. [p. 

The parameters appear in the n distribution amplitudes are listed below [3l|. We use the values at the scale 
fx = 1 GeV in our calculation. 



»2 m w 3 ru h 00 voo aw vw hoi h w 

0.25 0.015 -1.5 10 0.2 -3.33 -3.33 5.14 5.25 3.46 7.03 

It is reasonable to let M 2 = (3m p and M 2 = Pra^-^ , where j3 is a dimensionless scale parameter. Then we have 
u = m 2 /(m 2 + mlj and M 2 = I3m 2 m 2 j{m 2 + m^). 

From the requirement of the convergence of the OPE and the requirement that the pole contribution is larger 
than 40%, we get the working interval of the Borel parameter M 2 . The resulting sum rule is plotted with S01 = 
1.5,1.7,1.9 GeV 2 in Fig. [2] in the case of m 7I1 = 1.6 GeV. The sum rules for m 7ri = 1.8,2.0 GeV are similar. 
The numerical values of g p are presented here with their variations determined by the working interval of the Borel 
parameter 2.3 < M 2 < 2.7 GeV 2 and the range of the threshold 1.5 < s i < 1.9 GeV 2 . 



m 7Tl [GeV] 
ffpIGeV- 1 ] 



1.6 
2.5 - 3.2 



1.8 
2.6 - 3.3 



2.0 
2.6 - 3.4 



JW ! [GeV ! J 



FIG. 2: The sum rule for g p with m 7ri = 1.6 GeV, 2.3 < M 2 < 2.7 GeV 2 , and s i = 1.5, 1.7, 1.9 GeV 2 . 



III. SUM RULES FOR THE tti -> fm,bm COUPLING CONSTANTS 

Using the method outlined above we can further derive the sum rules for the ni —> fin, bin coupling constants. 
The effective Lagrangians of the process n\ — >• fin can be written as 

c h = -9/i x fw ■ Kg aP 

+gfd^n la x d a f ip ■ n + g} 2 d^n la x f w ■ d a n 

+gfn la x d a f ip ■ d n + gfn Xa x f 10 ■ d a d fi n . (15) 

Then the decay amplitude for this process is 

M {n x -> fx +n)= ig\ ( V ■ e*) + ig\ ( V ■ k)(e* ■ p) , (16) 

where g\ = g*~ — g 22 + g 2 ? — g£, e* is the polarization vector of the final fi meson, k and p are the momentum of 
the fi and the n\ meson, respectively. The interpolating current for the f\ meson is the axial-vector current 

Jfrfa) = ~y= [u(x)7 M 7 5 u(x) + J(x)7 M 7 5 d(x)] . (17) 

The fi meson couples to the above current through (0| J^ 1 (0)|/i(fc, A)) = ffxrrif^^. We consider the following 
correlation functions: 

U fl (k 2 ,p 2 ) = i J rfV^(,( 9 )|T{jA(,)J^(0)}|0) 

= -{9av - —o-jygpn — 



m 2 l * m 2 fi J (m 2 l - p 2 ){m 2 fi - k 2 ) 

[o},9^ + gjlk^Pu + gf^nPv + gfk^v + g}t%iA + 

, VctPv \ ( kpkf, ^ f^Jf l m h 

{g a u — ) {gpv 



m 2 i to 2 ^ {m 2 i — p 2 ){vn 2 i — k 2 ) 

x [g^g^ + gjlk^Pu + g}f(p^ - k^)p v + g}fk^(p v -k v ) + g}f(Pn - k^fo - k u )] h — 

= tow + A**] + < 18 > 

The axial- vector current (|17[) also couples to 7 = pseudoscalar meson 77/7/. However we can differentiate the 
contribution of the P-wave channel ni — > rj/r/'n to the correlation function from that of ni — > fin due to their 
different Lorentz structures. In a similar way to Section. [XXI we obtain the following sum rules for g 1 ^ and g 2 ± before 
the continuum subtraction: 

// 1 /. 1 ™/ 1 .9} 1 e fio ^ /M2+ " om - /M2 +--- 
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p Of 


Y 9 -A± 


2V2 I 


A da 3 


+ m 




\da 3 


da 1/ 



[ai] fdA±_ 8AA [a 



[cm] 

da2 J \ da 3 



da 3 da\ 



+ A±(u ,u , 0) - Aj_(u , u q , 0) + V|| (« , u Q , 0) + V|| (m , w , 



0) 



-V ± (u ,u ,0)-V ± (u ,u ,0) M 4 + K(uo) - ^K)]^ 2 ^ 2 + / B( u ) d M (,g 2 G 2 ) 
ffJ^m fl gle a ^ M2 ^<r/ M2 +--- 



(19) 



.An 



a[<*i,<Xi 
HI 



1 



l-Au 







ax]. 



2 



1 - 2u [ai] 
o I HI 



■ V 



(ff, 2 G 2 



(20) 



Here we omitted the terms ~ O(m^). After subtracting the continuum contribution, we get the sum rules for 

n ... . . r> . ■ ■■ ■ 



11U1U WU U1I11LLLX1 L11U LUlIlia ^^//i^y. i^LlLUl 11 LJ LI oXj Llllg Lilt: 

and g 2 fi , which are plotted with s i = 2.0, 2.2, 2.4 GeV 2 in Figs. HH 





s m =2.0GeV- 

s 01 =2,2GeV 3 

s 01 =2.4GeV 2 



Af-[GeVJ 



M 2 [GeV"| 



FIG. 3: The sum rule for gj i with 

m 7ri = 1.6 GeV, 1.7 < M 2 < 2.1 GeV 2 , and 

soi = 2.0,2.2,2.4 GeV 2 . 



FIG. 4: The sum rule for g 2 ^ with 

m 7ri = 1.6 GeV, 2.3 < M 2 < 2.7 GeV 2 , and 
soi = 2.0,2.2,2.4 GeV 2 . 



The derivation of the sum rules for g^ and g\ is almost the same as mentioned above. The interpolating current 
for the 61 meson reads 



4-> 



J 1 (x) = d(x) d ^u{x) 



(21) 



where d M = o> M — <9 M , and we define (0| J^ 1 (0)|6i(fe, A)) = fbi^u- There is another possible interpolating current for 
b\. J'Ax) — d(x)a t _ a/ u(x), which couples to 61 through (0| J' (0)|&i(fc, A)) = ifh l £^v P a^\k' T . The same current couples 
to the p meson through (0\J' llv (0)\p(k, X)} = ifj(^ki> — e*^). Unfortunately, the Lorentz structure of these two 
coupling will mix with each other in the correlation function 



Il'^(k 2 ,p 2 )=i I d*xe ik -*(ir(q)\T{J'^(x)j;^(0)}\0) 



(22) 



so that we are unable to separate the contribution of the b\ part from that of the p part. This is the consideration 
behind our choice of the interpolating current J^ 1 for the b\ meson instead of the tensor one. 

as 



The sum rules for g^ and g% read 



f b J,y bl e a °< /M2+Uor < /M2 +■■■ = 



216v / 2(ti u + rrid) 



b' a {u )-cj ) ' a {u )}M 2 (g 2 s G 2 ) 



(23) 
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U™1 



108v / 2(m u 



m d ) 



216w 



T(u ,u , 0) + T(uq, "o, 0) 



(dT 


dT \ laii 


(dT 


dT\ 


\da 3 


da 2 J 


\da 3 


da x ) 



M 



(1 - 2u )[<pcr(uo) + 4>a(u )] + UqUq^ (u ) - <p' a {uo)\ 



(24) 



They are plotted in Fig. [SJIH] after the continuum substraction. 




s l=2.0GeV- 

s ol =2.2GeV 2 

s ,=2.4GeV 2 



A/-[GeV-| 



FIG. 5: The sum rule for with m Wl = 1.6 GeV 
and soi = 2.0, 2.2, 2.4 GeV 2 . There is no stable 
working interval of M 2 for the sum rule. 



FIG. 6: The sum rule for g\ with m 7Tl =1.6 GeV 
and soi = 2.0, 2.2, 2.4 GeV 2 . There is no stable 
working interval of M 2 for the sum rule. 



The adopted values of the parameters ff 1 — 0.17 GeV and = 0.18 GeV 3 are obtained from Eq. (4.52) and Eq. 
(A. 20) in Ref. [32j, respectively. The extracted values of gj ± , g 2 x , g\ t , and g 2 x are collected in Table HI 



m ni [GeV] 


1.6 


1.8 


2.0 


M 2 


GeV 2 ] 


9/JGeV] 


-4.2 ~ -5.6 


-4.4 ~ -5.9 


-4.6 ~ -6.1 


1.7 


~ 2.1 


9% [GeV" 1 ] 


2.0 ~ 2.4 


2.1 ~ 2.5 


2.3 ~ 2.7 


2.3 


~ 2.7 


ffbJGeV] 


< 0.1 


<0.1 


<0.1 






9l [GeV- 1 ] 


1.9 


1.8 


1.6 







TABLE I: The numerical values of gj ± , g 2 ± , g\ x , and g 2 x . The working intervals of the Borel parameter M 2 are listed 
in the right column. Here "-" indicates the nonexistence of a stable working interval for the corresponding sum rule. 
The values presented in these cases are determined by M 2 = 2.0 GeV 2 . The range of the threshold is 

2.0 < soi < 2.4 GeV 2 . 



IV. SUM RULES FOR THE tti -> rpr,r]'n COUPLING CONSTANTS 

The interpolating currents for rj and 77' read as 

J"' = J" 8 sin (9 + J 711 cos6, 



where 



J^(x) 
J Vl (x) 



1 
1 

7! 



u(x)ij^u(x) + d(x)ij5d(x) — 2s(x)ij5s(x)] , 
u(x)i~/ 5 u(x) + d(x)ij 5 d(x) + s(x)ij 5 s(x)] , 



(25) 



(26) 
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and 9 = —19° is the mixing angle between ijs and the SU(3) singlet r\\. Their coupling to r\ and rjf are (0| J v (0)\rj(k)) 
X v and (0| J v (0)|r/'(fc)) = A ?) <, respectively. The correlation function involved in our calculation is 



U rl/Tl ,(k 2 ,p 2 )=i I d*xe ik -*(Tr(q)\T{Ji/n (x)J^(0)}\0) = <?G nM (k\p 2 ) + • • • , 
The sum rules derived for g n and g n ' are 



cos VUmt 



216v / 2(r7i u + md) 



f 2w\( —- — )+( — 



dT dT\ la21 

d^-d^) +T(uo,uo,0)+T(u o ,u Q ,0) 



M 



(27) 



(28) 



\ v 'f, l9v e u " r <' /M2+Uom ^ /M2 + ■■ 



sin f r,rTO; 



216v / 2(m u + m d ) \ [\da 3 da 2 
~ (f>' a {u )}(g 2 G 2 



[«i] 



dT dT 
da 3 dot\ 



+ T(u ,u , 0) + T(u , u , 0) 



M 



(29) 



We take ra, = 0.547 GeV, A^ = 0.23 GeV 2 , m v i — 0.958 GeV, and X v i — 0.33 GeV 2 in our numerical analysis [3£ 
The extracted values of g^ and g v i are presented here: 



m 7Tl [GeV] 


1.6 


1.8 


2.0 Af 2 [GcV 2 ] 




0.45 


0.45 


0.42 


ay 


-0.16 


-0.15 


-0.15 



PARTIAL DECAY WIDTHS 



It is straightforward to calculate the partial widths of n\ with the extracted coupling constants. The formulae for 
these partial widths are given in Eq. (|30]). 



9, 



96^771^ 

1 



P \ff |3 

y^r I ) 



r «?-**° 24W 



: ^ 192™^ 



i?;i 3 , 



(30) 



For completeness, we make rough estimates of the decay widths of some decay modes with the given coupling constants, 
although the corresponding coupling constants are extracted from sum rules without a stable working interval. We 
collected our results in Table UH together with the results obtained using other phenomenological models, e.g. Ref. 
fl9l | (PSS) and Ref. [lH (IKP). Here we also reproduce Table XIII of Ref. [35| to present the existing experimental 
results on the total decay width of k\ in Table Hill in order to compare with our predictions. 

The numerical values of the coupling constants of the modes tt\ — > pn, fiir are stable with the variation of m^. As 
a consequence, the partial widths of these two modes increase rapidly with m Tl due to their enlarged two-body phase 
spaces. 

Apparently, our results on the partial width of 7Ti are quite different from those obtained using Lattice QCD [40l ] 
and other phenomenological approaches such as the IKP and PSS flux tube model. In the flux tube model, the tt\ 
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m,^ [GeV] 




1.6 




1.8 




2.0 






IKP PSS this work 


IKP PSS this work 


IKP PSS 


Lattice 


this work 




[18] 


[19] 




[18] 


[19] 


[18] [19] 


[40] 




pn 


8 


9 


73 ~ 120 


12 


13 138 - 222 


16 




216 ~ 370 


fin 


14 


5 


69 ~ 122 


21 


9 96 ~ 175 


10.2 


90 ±60 


109 ~ 195 


blTT 


59 


24 


0.14 


62 


38 1.2 


43 


400 ± 120 


3.7 


r/n 








0.36 


.02 


.02 0.44 


.02 




0.45 


r/n 








0.02 





.01 0.02 


.01 




0.03 



TABLE II: The partial widths of the single-pion channels of tt\ in units of MeV, where IKP and PSS refer to the 
methods used in Ref. [l8| and Ref. [l9j . respectively. The partial widths for the channels bin and /i7r cited from 
IKP and PSS are simply the sum of their S'-wave and D-wave widths. 



Mode Mass (GeV) Width (GeV) Experiment Reference 



pn 


1.593 ± 0.08 


0.168 ± 0.020 


E852 


f5] 


7] 7T 


1.597 ±0.010 


0.340 ± 0.040 


E852 


f6] 


A* 


1.709 ±0.024 


0.403 ± 0.080 


E852 


\7] 


6l7T 


1.664 ± 0.008 


0.185 ±0.025 


E852 


m 


6l7T 


1.58 ±0.03 


0.30 ± 0.03 


VES 


[36] 


blTT 


1.61 ±0.02 


0.290 ± 0.03 


VES 


1 


b\TT 


~ 1.6 


~ 0.33 


VES 


[37] 


b\TT 


1.56 ±0.06 


0.34 ±0.06 


VES 


[38] 


/it 


1.64 ±0.03 


0.24 ±0.06 


VES 


[38] 


r\ 7r 


1.58 ±0.03 


0.30 ± 0.03 


VES 


[36] 




1.61 ±0.02 


0.290 ± 0.03 


VES 


[37] 


rj'n 


1.56 ±0.06 


0.34 ±0.06 


VES 


[38] 


blTT 


~ 1.6 


~ 0.23 


CBAR 


[9] 


pn 


1.660 ±0.010 


0.269 ± 0.021 


COMPASS 


fl] 


all 


1 fifi2+ 0015 


0.234 ± 0.050 


PDG 


[39] 



TABLE III: Reported masses and widths of the 7Ti(1600) from the E852 experiment, the VES experiment and the 
COMPASS experiment. The PDG average from 2008 is also reported. (Table reproduced from Ref. [Hj].) 



coupling to the two S-w&ve mesons is suppressed, leading to a small branch ratio of the mode ~k\ —> pn. One flux 

i Tl = 2.0 GeV is (in MeV) 

7T/0 : Tjn : tj'tt = 60 : 170 : 5 - 20 : - 10 : - 10 . (31) 



tube model prediction [4l| for widths for 7Ti with m - 

7r/i : 7r6i 



The quenched lattice QCD simulation also predicted a large width of the channel 6i7r [4(|, although the linear 
extrapolation approximation adopted there may lead to an overestimated width for this channel, as pointed out in 
Ref. [201 . However, the bin channel is severely suppressed in our calculation since the numerical value of g\ is found 
to be extremely small. As far as the channels r\n and rj'n are concerned, Chung, Klempt, and Korner argued that the 
channel n\ — > r\n is forbidden due to the requirement of Bose symmetry and J PC conservation in the limit that the 
j] is a pure SU(3) octet [42]. The tiny mixing between rj S and 771 should not reverse the widths of these two channels. 
In contrast, the width of r/n is at least one order of magnitude larger than that of the channel r/n in our calculation. 
Experimentally, the relative branching ratios for 7Ti(1600)'s three channels bin, rj'n, and pn are Q 

bin : rj'n : pn = 1 : 1 ± 0.3 : 1.5 ± 0.5. (32) 

In a summary on the VES results, Amelin [38[ obtained the relative branch ratios for the 7Ti(1600) as follows: 

hn : fin : pn : t/tt = 1.0 ± .3 : 1.1 ± .3 :< .3 : 1. (33) 
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VI. STRONG DECAYS OF THE ISOSCALAR 1 + HYBRID STATE 



Now we consider the strong decays of 7Ti, the isoscalar parter of m. We notice that I G J PC conservation restricts 
the possible decay channels to the S-wave tt 1 —> a 1 (1260)7r, fi(1285)ri, and the P-wave tti —> rjr]', 7r(1300)7r, 77(1295)77. 
The partial widths of the three P-wave channels are supposed to be relatively small due to their small phase spaces. 
In addition, the channel %% — > fiij is kinematically forbidden if the mass of tt\ is smaller than 1.83 GeV. Hence the 
dominant decay mode of 7Ti is tti — > a\n. 

The interpolating currents for the tti meson and the a\ meson are 



1 



1 

71 



uix^^uix) - d(x)j^ 5 d(x)] . 



(34) 



We define the coupling constants g^ and g\ x , similar to the coupling constants g^ and g 2 ± of the channel ~k\ — > fin. 
Apparently, the sum rules for and g% before the continuum subtraction are similar to Eq. (flUf and (|20p . 
respectively. The only difference between them lies in the hadron-level parameters, namely the masses and the 
overlapping amplitudes of the mesons involved in these two decay modes. Hence it is plausible to write down g\ x ss g 1 ^ 

and g\ x ps g 2 , if we simply adopt = f 7ri = 0.15 GeV 4 and f ai = f j t = 0.17 GeV in our numerical analysis. This 
leads to the estimate of the partial width of the mode fri — > a±n: T^ 1 ^ ai1T ps ST^ t -yf 11T . Here the factor 3 comes from 



the difference between the two channels' final states, namely n\ — > a\n ,a 1 n + , a^ir versus n® — > fin . 

If we denote the right side of the sum rules for g\ in Eq. (fT^j) and g 2 ± in Eq. (|2"0)) as R 1 ^ and R\ , respectively, 
then we have the following sum rules for g\ and g\„'- 



"V, ■■■), 
m n -> m a , ■ ■ ■ ) , 



(35) 



where we have ignored the SU(3) singlet octet mixing in the r\rj system and the factor e u o u om v /M because m 2 /M 2 < 
0.1 in our calculations. 







s 01 =2.0GeV- I 


s ,=2.2GeV 2 . 




s 01 =2.4GeV- 


1 2 3 


4 5 


M J [GeV 2 J 






- sqi =2.0GsV- 

- ioi=2.2GeV 2 

- i m =2.4GeV 2 



FIG. 7: The sum rule for gj, with 
m #1 = 2.0 GeV, 2.6 < M 2 < 3.0 GeV 2 , 
sqi = 2.0,2.2,2.4 GeV 2 . 



and 



FIG. 8: The sum rule for 
m #1 = 2.0 GeV, 2.6 < M 2 
sqi = 2.0,2.2,2.4 GeV 2 . 



9hv 



with 



< 3.0 GeV z , and 



We adopt /,, ( = 130 MeV and = 1.47 GeV [3l|. The mass of strange quark is taken to be m s — 0.15 GeV. The 
input parameters for the r\ light-cone distribution amplitudes involved in our calculation are as follows (jj, = 1 GeV) 



a,2 m w 3 m w 4 «oo «io vw h i h 10 

0.2 0.013 -3 0.5 0.2 -0.17 -0.17 0.17 0.26 0.15 0.38 
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The sum rules after the continuum substraction are plotted in Fig. [7][8l Their numerical values for 2.0 < Soi < 
2.4 GeV 2 are given here: 



[GeV] 


1.6 


1.8 


2.0 


Af 2 [GeV 2 ] 




-1.6 2.2 


-1.7- -2.3 


-1.7- -2.3 


2.6 - 3.0 




0.8 - 1.0 


0.9 - 1.1 


1.0- 1.2 


2.6 - 3.0 



Our calculation of 7Ti's widths are straightforward if we take advantage of the width formulae for n\ — > fiir: 
T ir«^a+ir- = 1 h*° (f/i ~* 9^, 9% -> 5ai . m h ~> rn ai ,m wi ->• m#J , 
r *?^/i»? = ^S-^/^r -> .9/irpff/i 9 , / 1 7 ) >w 7ri -> m^.m^ -> m,) . (36) 

The partial widths of 7Ti are listed in Table ITVl along with the results obtained using other approaches. As mentioned 
above, the only dominant channel for ffi is 7Ti — > ai7r whose width varies from 200 MeV to 600 MeV, depending 
heavily on the mass of %\ we adopt. 



m #1 [GeV] 


1.6 


1.8 


2.0 




IKP PSS this work 


IKP PSS this work 


IKP PSS this work 




[18] [19] 


[18] [19] 


[18] [19] 


aiir 


207 - 366 


72 28.2 288 - 525 


30.6 327 - 585 








8 10-19 



TABLE IV: The partial widths of tt\ in units of MeV, where "-" indicates that the corresponding channel is 

kinematically forbidden. 



VII. SEARCH FOR THE 1+ STATE AT BESIII 

Since the BESIII detector has an excellent photon resolution, it's very interesting to search for the 1 h state in 

m 2 — m 2 

the J/tp radiative process J/ip(ip') — > hi t o + 7. The photon spectrum peaks around E 7 = ^ h with a width 
rnhTh — (100 — 200) MeV. Such a process may be described by the following effective Lagrangian 

£ = coF^hlo + c' F^ a h h o^ (37) 

where F^ v is the electromagnetic field strength tensor, tp^ and h^ are the J/ip and 1 h field. Naively one expects 
the above branching ratio to be around 10~ 5 — 10 -4 . 

The isovector 1 h state tt\ can also be produced associated with other hadrons X at BESIII through the process 
J/t/j(tp') — > tti+X . For the production of the neutral component of 7Ti, the quantum numbers of X are I G = 1 + , C = — . 
Moreover, mx < mjM — rn- 7Tl — 1.5 GeV. From the above constraint, we get X = p°, b\, p(1450) if it is a single 
resonance or X = tt + it~ etc. Let's focus on the case X = p. Such a production may be described by the following 
effective Lagrangian 

C = cx^hl ■ fT + c 2 ^ht v ■ (T + ' ?™ + c^K a ■ % (38) 

Naively one expects the above branching ratio to be around 10 -4 — 10 ~ 3 . 

From Table mi the dominant decay modes of the isovector 1 h meson are pir, We urge our BESIII colleagues 
to search for m through the decay chain: J/ip(i//) — s- tti + 7, tt\ — >• pit —> n + ir~ir° or 7Ti — > /i(1285)7T°. /i(1285) is a 
narrow state with a width of 24.3 MeV. The /i7r° mode is also useful in the search of 7Ti, although /i(1285) mainly 
decays into multiple particle final states 4-7T, tjitit. The other important decay chain is J/tp(tf>') — > tt\+ p — > p + p + n — > 
2(7r + 7r _ )7r°. Once enough data are accumulated, one may also try to look for iri in the b±Tr, ijir, rj'w modes. 

The isoscalar 1 h state tti can also be produced associated with other hadrons X' at BESIII through the process 
J/ip(il)') — > 7i"i + X'. Now the quantum numbers of X' are I G = 0~,C = — . The possible candidates are X' = 
u), <j>, /ii(1170), w(1470), ir + ir~ir° etc. The ffi state mainly decays into aiir. Search of tti through the hadronic decay 
chain J/ip(ij}') — s- %% + u)/<fi —> ai + 7r + uj/<j> is challenging since it involves too many pions in the final states. BESIII 
collaboration may also search for fri through the radiative decay chain: J /^{ip') —¥ ffi + 7 — > ai + n + 7. 
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VIII. CONCLUSION 

We have studied the major strong decay modes of the J PC = 1 h hybrid mesons, including the isovector and the 
isoscalar cases. The coupling constants for these modes are extracted with the Light-cone QCD sum rule approach. 
Most of the sum rules obtained are stable with the variations of the Borel parameter M 2 and the continuum threshold 
sqx. For the other sum rules, we can not find a stable working interval of M 2 . 

Some possible sources of the errors in our calculation include the inherent inaccuracy of LCQSR: the omission of 
the higher twist terms in the OPE near the light-cone, the variation of the coupling constant with the continuum 
threshold soi and the Borel parameter M 2 in the working interval, the omission of the higher conformal partial waves 
in the light-cone distribution amplitudes of the pion or the 77, and the uncertainty in the parameters that appear in 
these light-cone distribution amplitudes. The uncertainty in the overlapping amplitudes between the interpolating 
currents and the corresponding final mesons is another source of errors. We merely give an estimated range for each 
involved coupling constant. The uncertainty of fh, which was not taken into consideration in this work, may broaden 
these ranges significantly. It's also understood that the a s correction may turn out to be quite large. We also omit 
the Oirn^) terms in the derivation of the sum rules for most of the coupling constants involved in our calculation. 

The partial widths of the 1 ^ hybrid calculated from the coupling constants extracted here are quite different from 
those obtained using other methods like the flux tube model and Lattice QCD etc. So far, the experimental data on 
the decay pattern of the 1 h hybrid is still not so accurate. We suggest the possible search of the isovector and the 
isoscalar 1 h hybrids in J/ip(ip') decay processes at BESIII. 
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Appendix A: The light-cone distribution amplitudes of the pion 

The 2-particle distribution amplitudes of the tt meson are defined as (3~i| 

(0\u(z) 7 ^ 5 d(-z)\ir-(P)) = if nPll f due l ^^(u)+ l -f w m 2 — Zt , f due^g^u) 

Jo 1 P z Jo 



(0\u(z)ij 5 d(-z)\Tr(P)) = f * m * C due lipz 4> p {u) 



m u + m d Jo 

(0\u(z)a afn5 d(-z)\iT(P)) = - l - U ™* ( Pa z p -p z a ) / due^ z ^{u) : (Al) 

3 m u + m d J Q 

where £ = 2u — 1, <f>^ is the leading twist-2 distribution amplitude, 4>( P ,<j) are of twist-3. All the above distribution 

amplitudes ip = {(j)^, (/> p , <j) a , g^} are normalized to unity: du <f>(u) = 1. 
There is one 3-particle distribution amplitudes of twist-3, defined as [3~fl ] 



{z)(J liylb g s G a f S {vz)d{-z)\TT (P)) = i {paP^p " PaPu9^ ~ PpP^dua + PftPv^t?) T(v,pz) , (A2) 



where we used the following notation for the integral defining the 3-particle distribution amplitude: 

T(v,pz)= I Vae-^ a ^ +va ^T(a d ,a u ,a g ) . (A3) 
Here a is the set of three momentum fractions a d , ot Ul and a g . The integration measure is 

J T>a = J da d da u da g 5(l — a u — a d — a g ) . 



(A4) 
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The 3-particle distribution amplitudes of twist-4 are 

(0\u(z)'f tl j 5 g s G a p(vz)d(-z)\n^(P)) = p^{p a zp - ppz a ) — f^mlA\\(v,pz) + (pfsg^ - Pa9p fJ ,)UmlA±(v,pz) , 

pz 

(0\u(z)'y f _ l ig s G a p(vz)d(-z)\ir-(P)) = p M (p a zp - p/3Z a ) — f„mlV\\(v,pz) + (p/sg^ - p a gf fi )f„mlV±(v,pz) , (A5) 

pz 



where G a p is the dual field G Q( g = ^e Q ,^ 7 5G 7 ' 5 . 

We also use the distribution amplitude given in Ref. [3lj ]: 



Mu) = 6u(l - u) (l + a 2 Cl /2 {0) , 
A(u) = 6uu 



1 + (1 + H a2 + 60% + y»fc)C£ /2 (0 + (-^02 - 6 m u 3 )Cl /2 (Z) 



(A6) 
(A7) 



16 24 „„ 20 
15 + 35 a2 + 207 ' 3 + T 7 ' 4 



1 1 

15 + 16 ~ 27 
18 



7 V3^3 -T^Vi) C% /Z (Q + 



27 



{-^o a2 -ik muj3 ) c * /2{0 



— r- a 2 + 21^4 {2u 3 (10 - 15u + 6u 2 ) \nu + 2u 3 (10 - 15u + 6u 2 ) lnu + uu{2 + 13uu)} , (A8) 



M(u) = g„(u) - </> n (u) , 

= 1 + (307,3 - %p\ ) c 2 i/2 (o + ( -3773^3 - ^ Pt - ^ p> 2 ) or (o 



27 2 81 2 



10 

,3/2, 



20 1 



,1/2, 



r (" ) = <i"i 1 - " ) { J + j •"'//. t - ~ 773^3 - - ~ p*a 2 1 r j "!4S 

1 



T(a) = 360?73Q;„ada 9 |l + w 3 - (7a g - 3) J- , 

V||(a) = 120a u a d a g (u O o + v w (3a g - 1)), 
= ^0a u a d a g a w (a d ~ a u ) , 



Vi(a) = -30a 



/i 00 (l - a„) + /i i a g (l - a g ) - 6a u a d 



'10 



"4-L(a.) = 30a 2 (a n -a d ) 



/ioo + ^oia 9 + - h 1Q (5a g - 3) 



(A9) 
(A10) 

(All) 

(A12) 

(A13) 
(A14) 

(A15) 



where C™(£) are Gegenbauer polynomials. 

The definitions and the specific forms of the rj light-cone distribution amplitudes adopted in the text are similar to 
those of the pion. For more details see Ref. |3l| . 



